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Abstract. We review the appearance of multiple scalar fields in linearized SFT based cos- 
mological models with a single non-local scalar field. Some of these local fields are canonical 
real scalar fields and some are complex fields with unusual coupling. These systems only 
admit numerical or approximate analysis. We introduce a modified potential for multiple 
scalar fields that makes the system exactly solvable in the cosmological context of Friedmann 
equations and at the same time preserves the asymptotic behavior expected from SFT. The 
main part of the paper consists of the analysis of inhomogeneous cosmological perturbations 
in this system. We show numerically that perturbations corresponding to the new type of 
complex fields always vanish. As an example of application of this model we consider an ex- 
plicit construction of the phantom divide crossing and prove the perturbative stability of this 
process at the linear order. The issue of ghosts and ways to resolve it are briefly discussed. 
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1 Introduction 

A new promising class of cosmological models originated from string field theory (SFT) [1] 
and p-adic string theory [2] has emerged recently and attracts nowadays a lot of attention 
[3-9]. The main distinguishing feature of such models is the presence of specific non-local 
operators in the action. There is a strong theoretical reason to consider such non-local 
interactions. In order to have a good effective model the UV-completeness of the underlying 
theory is essential and these non-local operators are crucial for making SFT and p-adic string 
theory UV-complete [1, 2]. Furthermore, such non-local operators have proven to give rise to 
models with new appealing properties in the context of cosmology, 1 and various interesting 
results have been obtained recently. 

In application to the Dark Energy (DE) problem, one can easily construct models obey- 
ing the phantom divide crossing. This is of interest for the present cosmology. Indeed recent 
results of WMAP [11] together with the data on la supernovae and galaxy clusters measure- 
ments, give the following bounds for the DE state parameter wde = — l-02+Q';[g. The crossing 
of the phantom divide as well as the phase with w < — 1 are not excluded experimentally. 
Moreover, the experimental data for cosmological epoch corresponding to redshift z > 1 leave 
even more freedom for possible values of w during that phase. 

The phantom divide line vu = — 1, corresponding to pure vacuum energy, gives the 
lower bound of the Null Energy Condition (NEC), whose physical motivation is to prevent 
instability of the vacuum. Therefore in order to obtain a stable model with w < — 1 one should 
construct an effective theory, with the NEC violation, from a fundamental theory which is 
stable and admits quantization. This is a hint towards SFT inspired cosmological models. 
It is therefore important to find models where the DE crosses the phantom divide. Among 
other cosmological models with w < — 1 and free of instabilities that have been constructed 
we can mention the Lorentz- violating Dark Energy model [12], the ghost condensation model 
[13], model with the kinetic gravitational brading [14], model having the dust of the Dark 

1 Other non-local models were considered as well in the literature [10] in the cosmological context. 
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Energy [15], NEC violating models based on the Galileon theories [16], and the brane- world 
models [17]. 2 

Other results from non-local cosmological models have been obtained in the context of 
inflation [4, 5]. Such models have in general the remarkable property that slow roll inflation 
can proceed even in presence of an extremely steep potential [6]. Furthermore they can 
produce non-gaussian signature for the cosmic microwave background. It has been shown in 
[4] that the non-linearity parameter which characterizes this non-gaussianity can be large, 
as opposite to what is described by the canonical single scalar field inflationary models, and 
observationally distinguishable, for instance, from Dirac-Born-Infeld inflation models. 

Recently an attempt to build a non-local generalization of Einstein gravity has also 
been developed in the literature [7] (see also [21] for an application of the diffusion equation 
method to non-local gravity models). One result obtained in this direction is the construction 
of an exact solution which describes a non-singular bounce at the origin and a de Sitter phase 
in the large time asymptotics. 

The main focus of this paper is on a particular class of cosmological models that nat- 
urally arise in certain regimes of SFT [3-9]. As it will be explained in details in section 2, 
such models effectively contain (infinitely) many scalar fields, some of which are characterized 
by the unusual feature of having complex masses squared. Solutions to this type of models 
have been analyzed from the cosmological point of view and have proved to be interesting in 
order to address the questions of the non-gaussianity and the phantom divide crossing. It is 
however vital to consider the evolution of the (at least linear) perturbations. This step has 
not been yet covered in full generality for inhomogeneous perturbations in the presence of 
fields with complex masses squared. It has been explored only to some extent in [6, 8, 22, 23], 
where space homogeneous perturbations were analyzed. The goal of this paper is to fill this 



The paper is organized as follows. In section 2 we review the general properties of this 
class of SFT derived models. In section 3 we present the specific model that we are going 
to consider and we show that it admits a class of simple exact solutions in the cosmological 
background represented by spatially flat FRW universe. The general equations describing 
the evolution of cosmological perturbations around these exact background solutions are 
presented in section 4. In section 5 we specialize to the case corresponding to a truncation 
of the local system to a system with two complex conjugate fields with complex conjugate 
masses. In section 6 we undertake the numerical study of the perturbations and describe the 
results obtained. Section 7 describes the issue of ghosts in our model. 

2 SFT based non-local models 

Non-local models arise naturally in SFT [1] where, under certain approximations, one gets 
the following space-time action 



2 A remark on stability is necessary. Ghost condensate models look unstable and this is cured by adding 
higher derivative terms. This in turn implies the presence of other ghosts and related instabilities [18, 19]. 
Nevertheless, since a cutoff is assumed and the wavelength of this ghost is smaller than the cutoff, the 
instabilities can be considered irrelevant. The stability claim therefore relies on the fact that a well behaved 
UV-completion is assumed to exist. This is actually the case in [12], but it seems that ghost condensate 
models cannot be UV-completed in a usual way [20]. 



gap. 
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Here J-q{z) is an analytic function on the complex plane with real coefficients in the Taylor 
expansion. $ is a scalar field from the string spectrum and can be for instance the string 
tachyon. To be more specific consider for example cubic SFT, which may be either bosonic 
or fermionic. Performing standard SFT calculations one arrives to the action for space-time 
fields. This consists of two parts: the quadratic part, which is the kinetic term, and the cubic 
interaction among all the fields. Schematically we can write it as 



where f3 is a parameter determined exclusively by the conformal field theory and □ is the 
d'Alambertian operator. It is natural to expect /3 < 0, corresponding to convergent Feynman 
graphs at large momenta. From this specific example we see the non-local character of the 
action that is a general feature of SFT based models. 3 If furthermore one integrates out all 
the modes of action (2.2) except one and if only the low mass excitations of the SFT are 
retained, one obtains 



as an effective action. With a simple field redefinition $ = e" eft one obviously gets an 
action of the type outlined above. 

The presence of the operator J-q is the key modification with respect to a canonical 
scalar field action and its properties would play the most important role for the resulting 
model. In SFT this operator can be computed perturbatively by means of the so-called "level 
truncation scheme" [1]. In general the action of the form (2.1) is just an approximation. In 
principle a more involved entanglement of non-local form factors and scalar fields would be 
obtained. An example of the operator J-q is given by what one obtains form the tachyon 
field in the low-level approximation, where one may have J-q = (□ — m 2 )e 2l3n . A specific 
expression for JFq does not however simplify the following analysis, especially in the case of 
a curved background, and therefore we keep it general. Moreover, this operator may get 
modified if we couple our model with other ingredients. Since this may result in possible 
restrictions on the operator and coupling structures in the theory, it is important to keep J-q 
general. 

From the SFT point of view, the full potential V = — ^^(O)^ 2 + Vi n t in (2.1) is expected 
to have at least two distinct extrema. The scalar field exhibits a non-perturbative transition 
from one extremum to the other and eventually stops. The best understood processes of 
this kind in the SFT dynamics is the so called rolling tachyon (see, for instance, [24] and 
refs. therein), described as a time-dependent kink-like solution. Several solutions in flat 
background were constructed numerically [24] and in some approximations analytically [25] . 
Taking the linearization of the potential around a given extremum, say <!> = <l?o + V 9 ; a 
significant analytic progress can be achieved. The linearized action becomes quadratic in the 
scalar field ip 



In the linear approximation, the action will actually have such a form and no other structures 
are expected. 

3 Appearance of higher derivatives is not an exclusive feature of this theory. Non-commutative theories, for 
instance, also have higher derivatives, but these non-local structures are very different. 




(2.2) 




(2.3) 




(2.4) 
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For the action (2.4) one can establish an equivalence of Lagrangians [9] (notice that 
this mapping is possible only for the quadratic linearized action and not in presence of a 
potential, like in (2.1)) 



where in the r.h.s. the Jj are roots of the characteristic equation J- (J) = and there are as 
many terms in the series as roots. In other words there is a mapping of a non-local theory 
with one non-local scalar field to a local theory with many local scalar fields Tj. To each root 
is associated one local field, thus in the case of many but finite roots one obtains a genuine 
local theory with many scalar fields. When infinite roots are present, the non-locality of 
(2.4) is rephrased in the appearance of an infinite number of local fields. In general the 
roots Ji may be complex and, in this case, they appear in complex conjugate pairs. As a 
matter of fact this is not a problem for the theory, since the local fields are non-physical. 
They only provide a different mathematical description. It is the original non-local field 
and the related quantities that must be real, since they represent physical excitations. This 
translates in requiring that the hermiticity of the local action is preserved when a pair of 
complex conjugate roots is present. This means that the fields Tj also enter the action in 
complex conjugate pairs. The point we want to emphasize is that instead of the usual rr* 
quadratic combination we have terms of the form 



This quadratic expression in general cannot be diagonalized in terms of real fields and this is 
one of the distinguishing features of such models compared to canonical ones. It is however 
simpler to work with (infinitely) many local fields rather than one non-local field. 

The only thing which is important is the particular vacuum around which the action 
is linearized and the associated spectrum of the theory. There are many different possible 
vacua in string theory and the class of models we deal with practically distinguishes these 
vacua through the form of the operator T . The latter is in turn characterized by its roots Ji 
and there may be vacua in which we have no zeros at all (it is believed that the true tachyon 
vacuum is of this kind with, for example, J- = e 2 ^ n ). 

3 Exactly solvable model 

Looking for a cosmological scenario and expecting a similar behaviour in a slightly curved 
background, one is lead to consider a scalar field that interpolates in between of distinct 
vacua of the potential in a curved background and eventually stops (i.e. ip tends to 0). Then, 
in the vicinity of the final vacuum, one can construct approximate background solutions 
for a model where many scalar fields are minimally coupled to gravity and all Tj tend to 
[9]. Furthermore, it was explicitly shown in [22] that exists a simple generalization of 
such a model with (infinitely) many scalar fields minimally coupled to gravity such that the 
equations of motion in the spatially flat Friedmann-Robertson- Walker (FRW) metric become 
exactly solvable. This generalization consists of an additional quartic potential 



i 



niu-j^n + rtin-r^r*. 
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that makes the full equations exactly solvable for specific values of the parameters on. (It was 
argued in [22] that such a modification is obviously the next to the quadratic order expansion 
of the original potential around an extremum in the cases where the potential is even around 
this extremum, i.e. in cases where the expansion starts with V = \m 2 tp 2 + jyA^ 4 + ...). Such 
an additional term vanishes more rapidly than the other terms of the original model when 
evaluated on solutions with the scalar fields going to zero. Therefore, in the far asymptotic 
regime, the results for the model with and without the quartic extra term should coincide. 

To be precise, the model of primary concern in the present paper is described by the 
following action 



16ttG 



N 



55 (- £ \ («r<w^ 



where in addition to the quadratic piece a specific potential of the fourth degree is present. 
We work in 1 + 3 dimensions with the signature (— , +, +, +), the coordinates are denoted by 
Greek indexes fi, v, . . . running from to 3. Spatial indexes are a,b, . . . and they run from 
1 to 3. Gn is the Newtonian constant, 8irG]y = 1/Mp, where Mp is the Planck mass, g Q 
is the open string coupling constant, G = G^ jg 2 . 9uv 1S the metric tensor, R is the scalar 
curvature, A is a constant. The fields are dimensionless, while [g ] = length. 

We impose a specific relation between the parameters, namely Jj = — ai(pn + 3Hq), with 
H = y/ 87rGjvA ^ guch a relationj 

even not being expected to hold always, is the one that 
one gets in SFT considering the tachyon vacuum in the low-level approximation and with 
dS background metric [1, 26]. Originally, because of the connection with SFT models, one 
would set the Jj as the fundamental family of parameters. On the other hand the modification 
introduced here emphasizes the role of on as opposite to Jj. Consider for example one specific 
Jj. In general there are two possible a associated to it. If one wishes to include both of them 
in the action one realizes that, due to the specific form of the quartic potential, the two fields 
corresponding to the two different a are really independent. For instance they cannot be 
introduced through their linear combination, as opposite to the case of free fields with the 
same mass. They simply need to be introduced in the action as two general fields, with the 
special feature that the associated J turns out to be the same. Moreover, as it will become 
clear, the exact solution is really parametrized by a. Finally, as it will be discussed, for our 
purposes we are specifically interested in a particular range of values for a. This practically 
associate in an automatic way one single oti to each Jj in most of the cases. 

In the sequel we are going to explore the case of spatially flat FRW cosmology, with our 
scalar fields system minimally coupled to the background metric 

g^ = diag(-l, a 2 (i),a 2 (i),a 2 (i)) , 

where a(t) is the scale factor, t is the cosmic time and we shall use the Hubble parameter 
H(t) = a(t)/a(t), denoting with the dot the derivatives w.r.t. the cosmic time. In such a 
background the modified local action has the following equations of motion 

3H 2 = A-kG ( f i + J i T i) + 37rG (^2 aiT ^j J + 87tG ^ A » 



(3.2) 



- 5 - 



and 

fj + 3Hii + JjTj + 67rGajTj a j r j — > f° r au * • (3-3) 

j 

It is easy to check that they admit the following exact solution 



n = T i0 e ait 



H = H - 2vrG ^ a^e 



2a t t (3.4) 



This solution is valid for any number of fields (including a single field and infinitely many 
fields) and for any values of the parameters 014 (i.e. real or complex). We see that if Re (a j) < 
for all i then solutions vanish and moreover the quartic term in the potential practically 
vanishes at sufficiently shorter times and we are left with free fields. Thus for large times 
the model without a quartic potential is restored and we return to more general SFT based 
models without requiring any other conditions on the potential to be satisfied. More generally, 
(3.4) is a solution of linearized Friedmann equations in any non-singular potential with 7$ = 
0, H = Hq as a stable fixed point, according to Lyapunov theorem, if Re{oti) < 0. In our 
particular setting we deal with a potential for which this solution is exact for the full action. 
This would simplify the task of considering perturbations, but a generalization of the results 
to other forms of the interaction seems straightforward. 

From the physical point of view, the solutions above describe a very interesting phe- 
nomenon, the phantom divide crossing. If there are complex roots Jj, H oscillates around 
the values Hq and the equation of state parameter w = p/p crosses the value w = — 1 (p and 
p are the collective pressure and energy density). Indeed 

2 H 

and H is an alternating function. Notice that such a crossing is generically not possible with 
single ordinary scalar field. With two scalar fields one can obtain such a behavior, but one of 
the two fields must be a phantom (ghost) in this case. In our setup, in presence of complex 
roots, we also have at least two fields, which are however just an effective description of the 
original theory with a single non-local field. Moreover, recall that the complex fields enter 
the quadratic part of the Lagrangian in the following way: 

nip- j^n + rtip- J*)r*. 
This can be recast in terms of real components as 

XPXi ~ ipiOtpi - mi(xi ~ + ^niXi^i , (3.6) 

where r, = ^(Xi an( i Ji = nii + irii. This Lagrangian would be a partial case of the so 

called B-inflation model [27] provided the "mass" rrii becomes zero for one of the fields, say 
ipi. Also, this resembles the so-called quintom model [28] where two fields (one normal and 
one phantom) are used as well, but without a quadratic coupling. Because of this coupling 
we cannot exactly identify the states as normal or tachyons or ghosts, since these notions 
refer to diagonalized Lagrangians with real fields (in our case we cannot make it diagonal 
with real fields). This does not mean however that the issue of ghosts is not present in our 
model. We will return to this question in section 7. 
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4 Cosmological perturbations 



We want to consider the problem of cosmological perturbations for the exact solutions pre- 
sented in the previous section. We focus on the scalar type perturbations. These are induced 
by energy density inhomogeneities and are the one that may cause cosmological instabil- 
ity, as opposite to vector type perturbations, that are known to decay rapidly and tensor 
type perturbations, that correspond to gravitational-wave perturbations of the metric and 
do not induce any perturbation in the perfect fluid at the linear order. The three kind of 
perturbations are decoupled and thus one can safely consider them separately [29]. 

We start presenting the equations for spatially inhomogeneous perturbations of the 
modified action introduced in the previous section. We refer the reader to [29-31] for the 
details on the general derivation and to [8] for the conventions adopted here. In our setup the 
primary quantities are the scalar fields Tj. It is therefore more transparent to use the fields 
themselves rather than the notions of energy and pressure for the corresponding components 
of the stress-energy tensor. It is also important that the spatial dependence of the perturba- 
tion variables can be represented as e tkaXa for the spatially flat FRW universe, where k a is a 
3-dimensional vector and k 2 = 5 ab k a kb- With this representation all the relevant information 
about the space-inhomogeneities are collected into the comoving wavenumber k. 

After a number of tedious manipulations the relevant equations for scalar perturbations 
read [8, 31] 



Ti T<i\ ( „ ,', k 2 

a- 



)- 


-( 




Ti \ 



T, 2 



— 7- yjj + GirGaj ^2 a i T i 
+ UttGS^ a k T k f k ( a^Cik ~ aj^-Cjk ) , 




i + eH(2-6w + 3cf) + e[H(l- 3w) - IhwH 1 + 9H z c z s + 
k 2 l 2 



0,2 pp + p t7 



^4 + 6nGa k ^ ayr 2 T k f k ffC,ki 



where Cij = — is the gauge invariant variable for the scalar fields perturbations and 
e is the gauge invariant total energy density perturbation. In writing the equations we have 
introduced the definition c 2 s = p/p for the speed of sound. 

Equipped with these general equations we can consider as a background the exact 
class of solutions (3.4). Substituting the expression for Tj, one gets some cancellations and 
simplifications. Keeping track of all of these, the final result is: 

m ! / x k 2 (ol — a-:) 



a 1 H 

\ (4-1) 

rWkZ 2akt {Cik + Cjk) + 3H 2 e , 
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/ H h\ ( H 2 H k 2 ^ 

\ H HI \ H 2 H a 2 



k 2 4 (AirG) 2 \ - ^ „.2„.2_2 _2 J2a k t„2a,t 

a 2 3 HH 2 

k,l 



(4.2) 



where the equations of state parameter and the speed of sound have been expressed in terms 
of the Hubble parameter and its derivatives, according to Friedmann equations. 

A general analysis of the system given by (4.1) and (4.2) is difficult, even in the case 
where the background solutions have the simple form (3.4). Many models with multi scalar 
fields derived from non-local model with a single scalar field and the related problem of 
cosmological perturbations have already been considered in the literature. Nevertheless, for 
all the motivations discussed above, the model presented here represent to a certain extent 
a novelty in this context. In the analysis of the cosmological perturbation we will therefore 
concentrate on the simplest case that capture the interesting features of the model presented 
in this work, namely the case of two complex scalar fields. 

5 Two complex conjugate roots 

When considering complex solutions one should remember that the physical quantities that 
appear in the original model must be real. In the case of two complex conjugate solutions, 
labelled by the parameters a\ and «2 = ck*, that we wish to consider, this amount to choose 
the integration constants in such a way that t\ = r| . This is exactly the setting we are going 
to consider. Using the representation a\ = a\ = a = — | + i\ or a = |a|e"^ a and using 
the explicit solutions for t\ = t% = Toe * one can write all the quantities appearing in the 
differential equations by means of trigonometric functions. Upon substitution and after some 
algebra the equations of interest take the form: 

/ £2 \ yeH 2 

X + X (3-ff + x + y tan(yi)) + x 3x# + yx tan(yt) + — = — - , (5.1) 

\ or / cos(yr) 



e + i \ + 4^ + x + y tan(yt) j 

/ if- 1.2 \ l2 

+ e[ 6H 2 + UH + 2—f + 3xH + 3yHt&n(yt) + ^ = - TT9 = , - . (5.2) 
\ H z a z I H z a z cosyyt) 

Note that all the quantities appearing in the equations are real. All the arbitrary phases 
have been eliminated. This can be done in all generality corresponding just to a shift in 
the time variable. The function x is a ls° rea l an d is related to the original perturbation 
parameter ("12 for the scalar fields as x = i 8nG ^ a 3 ^ T °^ e~ xt Ci2- Notice that one can always 
take y > without any loss of generality. In fact one sees from the explicit form of the 
equations that a change of sing for y can be always reabsorbed as an extra minus factor 
in the definition of x- This is not important if we are interested only in separating regions 
in the parameters space where the perturbations are confined from regions where they are 
not. Moreover, y just sets the periodicity of the functions appearing in the differential 
equations and can be therefore conveniently rescaled when studying the numerical evolution 
of the perturbed system. This corresponds to a change of variable from the cosmic time t 
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to the dimensionless time coordinate T = yt. With such a rescaling of the time variable 
the perturbative equations are characterized by a set of dimensionless parameters, Ho/y, tq, 
x/y and k/y, once an extra 1/y factor is reabsorbed in the definition of x- Note that with 
this redefinition of the time variable one also has H{t) = yH(T). From a practical point 
of view the rescaling is completely equivalent to set y = 1 in (5.1) and (5.2), with a set of 
parameters, Hq, tq, x, k, that are now dimensionless. We therefore adopt this choice, setting 
y = 1 hereafter and we still denote the time variable as t. The explicit expression for H in 
the parameterization adopted is: 

H = H + h = H - 47rG\a\\T \ 2 e- xt cos(t - <j) a ). (5.3) 

This makes transparent that for certain parameters one can have an expanding universe 
regime, H > 0, with H reaching the constant value Hq at large times. It is also useful to 
keep in mind that 

H = -8^G|a| 2 |ro| 2 e" xi cos(t) 

and 

a = a G exp (H t - 2irG\T \ 2 e~ xt cos(t - 2<f> a )) . 

We stress one more time that hereafter all the quantities appearing in the equations and in 
the definition of H and a must be understood as rescaled dimensionless quantities. Note 
also that we emphasized that we are particularly interested in situations where the quartic 
potential vanishes for large times, corresponding to x > 0. As discussed, the rescaling does 
not arm this point and we can safely work in terms of positive rescaled dimensionless x. 

The system of equations (5.1)-(5.2) is in principle ready to be analyzed. An analytic 
solution for this problem is however not easy to find. In order to proceed with a numerical 
analysis it is worth noting that there are some special points in the evolution of the system 
under consideration. Namely, when t = ir/2 + mr one has divergent coefficients appearing 
in the differential equations. Indeed at these points cos(£) becomes zero. The numerical 
integration process may encounter some problem when it reaches these points if the precision 
is not high enough. Even worse, the system could be really singular. Expanding in series 
around such points the differential equations read: 

tX-X~x X = -(-l) n Hle, 

k 2 

ts-e- 3H n e = {-) n jpr^X , 

where the constants H n and a n are respectively the functions H{t) and a(t) evaluated at 
the n-th singular point. From (5.4) it is quite clear that nothing happens close to the 
singularity. Numerical integration of this equations and of the non-approximated ones close 
to the singularity shows that one can pass the singular points without any problem and, if 
the precision is high enough, all the features of the differential system are captured. 

6 Numerical study 

We describe here the numerical results obtained starting form the system of equations pre- 
sented in the previous section and the picture that emerges for the evolution of the pertur- 
bations when the different quantities characterizing the background and the perturbations 
themselves are varied. We particularly focus on distinguishing between situations where the 
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(5.4) 



system remains stable and situations where the perturbations grow making the system un- 
stable. For simplicity we set the constants G = l/8ir and ao = 1. This choice does not affect 
the evolution of the system in any way. 

In order to numerically integrate the system of differential equations (5.1)-(5.2) exact 
background solutions must be specified. In the case of two complex conjugate roots consid- 
ered here, this corresponds to specify three real dimensionless parameters, Hq, tq, x. We 
specifically focus on the cases of asymptotically vanishing quartic potential, corresponding to 
positive x parameters. Within this class of solutions, the Hubble parameter H evolves with 
damped oscillations around the constant value Hq. In the analysis of the perturbed system 
we concentrate on the case of cosmologies with expanding universe. This corresponds to ask 
H > in all the stages of the evolution. Such a request, in the region where Re(a) < 0, 
is easily fulfilled with a balanced choice of a and tq with respect to Hq. In figure 1 is 
shown the typical evolution of the Hubble parameter considered in our study of cosmological 
perturbations. 




Figure 1. In the left graph, the typical evolution of the Hubble parameter considered in our setting is 
plotted against the e-folding number. It evolves with damped oscillations towards the asymptotic value 
H , never crossing the zero. In the right graph, the corresponding evolution for the equation of state 
parameter w is shown. The crossing of w = — 1 that characterizes the case of complex roots is evident. 
The value of the constant parameters considered in the plots are: Hq = 0.07, x = 0.01, tq = 0.25. 

One comment is necessary in this regard. The constant of integration To can be in 
principle tuned at will for the purpose of reproducing the desired behaviour for H. However, 
from a physical point of view, this choice is not completely arbitrary. Equation (3.2) shows 
indeed that the value of tq is an important parameter in determining the origin of the Dark 
Energy content of the model. 

One additional parameter that characterizes the perturbed system and that must be 
specified in order to numerically study the evolution of the perturbations is the comoving 
wavenumber k. This always enters the differential equations in a specific way. Namely, it 
always appears in the combination k 2 /a 2 and it is therefore mainly relevant in the early stages 
of the evolution. At large times it become less and less important, due to the suppression given 
by the increasing scale factor. In particular this combination enters the non- homogeneous 
term of the differential equations and is therefore an important parameter in setting the 
coupling between the scalar fields perturbation parameter and the energy density one. Note 
however that it appears in a non-symmetric way in the differential system (5.1)-(5.2). It only 
appears in the non-homogeneous term of (5.1), while it enters in a similar way into the left 
hand sides of the two differential equations. 

Finally a set of initial conditions for the differential problem must be specified in order 
to numerically integrate the differential equations. In particular one would like turn on a 
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small initial non-zero value for one or both the quantities e and x> a t some given time. The 
initial value problem requires also a set of initial conditions for the first time derivative of 
the two perturbative parameters. 

There are mainly two different situations emerging from the study of the differential 
system (5.1)-(5.2). In both cases one observes that the perturbations decrease to zero at 
sufficiently large times, giving in principle a well behaved evolution for the system, such 
as the one shown in figure 2. There are however particular regimes where one must pay 
attention. While in general one could say that the perturbations evolve with oscillatory 

x 




Figure 2. The evolution of e and \ f° r Ho = 0.07, x — 0.01, tq = 0.25, k = 0.5 is shown. The initial 
conditions are: e(0) = 10~ 4 , e(0) = 0, x(0) = 10~ 3 and x(0) = 10" 4 . 

behaviour according to various possible different paths, but eventually always decreasing to 
zero value for both e and x> i n some particular cases one observes a fast initial growth for e. 
Precisely, one does not observe divergence for e. On the contrary, it eventually vanishes, but 
in the initial stages of the numerical evolution it can exit the perturbative regime. 

Indeed when one has non-zero initial conditions for x> irrespectively of the initial value 
of e and its derivative, one naively observes that if Hq < 1 and if k 2 /Hq is of the same order of 
magnitude of 1/x, for values of x taken in the initial stages of the evolution, e increases. If the 
non-homogeneous term in (5.2) is large enough, the energy density perturbation parameter 
will exceed the range of values where the perturbative regime makes sense. An example of 
such a situation is depicted in figure 3. This means that if we start with non-zero initial 




Figure 3. The evolution of e and \ for Hq = 0.07, x = 0.01, tq = 0.25, k — 15 is shown. The initial 
conditions are: e(0) = 10~ 4 , e(0) = 0, x(0) = 10~ 3 and x(0) = 10~ 4 . Even though both x and e stay 
confined and eventually decrease to zero value, e goes trough an initial phase where it shows a fast 
growth and exits the range of validity of the perturbative regime. 

conditions, one can always find a value of k large enough to trigger a fast and large growth 
for e. In fact the conditions naively identified for such a phenomenon to occur do just mean 
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that the non-homogeneous term should remain significantly large in the initial stages of the 
evolution. The scale factor appearing in the non-homogeneous term, on the other hand gives 
an exponential suppression of the latter with a rapidity roughly set by Hq. This can be 
always compensated with an appropriately large value of k, so as for the balance of the other 
factors appearing in the non-homogeneous term of (5.2). The exponential suppression of the 
scalar field solutions, due to the positive values of x considered here, also play a role in this 
subtle interplay of k with the other parameters. In fact the amplitude of the oscillations 
covered by e increases roughly in a polynomial way with k/x. Specifically, for small values 
of k/x it increases like (k/x) 2 . When the ratio k/x gets larger this power law behaviour is 
smoothed out and the maximum range of values spanned by e increase more slowly with k/x. 

In the opposite case, where one sets zero initial conditions for x an d X, the spectrum of 
perturbations goes to zero without any particular surprise, s evolves with different oscillatory 
modes, spanning a range of possible values of the same order of magnitude of the initial 
conditions and eventually vanishes. In these cases x is a ls° excited, but only up to values 
that are practically negligible. Such an asymmetric behaviour dictated by different initial 
conditions shows that exists a strong transfer from x to e, but not in the opposite direction. 

One can therefore conclude that for any value of the parameters Hq, tq, x and for 
any initial conditions, provided that x an d X have not simultaneously zero initial value, 
it is always possible to find a k large enough to obtain a rapid and large growth of e. 
This can bring e out of the range of validity of the linear approximation adopted in the 
study of the cosmological perturbations. When this happens one cannot relay anymore 
on the perturbative computations. Since, however, the perturbations remain confined and 
eventually they all go to zero, it is an interesting question whether this phenomenon signals 
a real instability of the system or if a more careful analysis, such as considering higher order 
perturbations, could reveal that the system is indeed stable irrespectively of the value of k. 

As it was explained at the end of section 3, the case of (two) complex conjugate roots 
corresponds to a situation where the phantom divide crossing is realized. The analysis of this 
section shows that for a large range of parameters the system has a stable behavior under 
inhomogeneous perturbations and the perturbations are in any case confined and vanish. 
This allows us to claim that the exact solutions presented in this paper give a stable (at least 
up to linear perturbations) model with the equation of state parameter oscillating around 
the value w = — 1. 

7 Taming the ghosts 

The numerical analysis of the previous section shows that at the linear order in perturbations 
the system is stable for a wide range of parameters and initial conditions. Even in cases where 
the perturbations grow, posing a question on the validity of the perturbative approach, it 
is enough to wait a sufficiently long time to see that perturbations finally go to zero. One 
could therefore conclude that there are no divergent modes in perturbations. This however 
does not mean that ghosts, which are effectively present in our setup, would not trigger 
instabilities at higher orders in perturbations or when coupled to other fields [32-34]. It is 
therefore important to understand to what extent ghosts do not pose a threat for the system 
and our results are applicable. 

To clarify the picture we return to the SFT motivation of the model we analyzed. 
We have considered the next-to-the-linear approximation of SFT dynamics around a given 
vacuum. There are many different possible vacua in string theory and the class of models 
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we deal with practically distinguishes these vacua through the form of the operator T . The 
latter is in turn characterized by its roots Jj. There may be vacua in which we have no 
zeros at all (it is believed that the true tachyon vacuum is of this kind for example). Such a 
vacuum is definitely ghost-free. There may also be vacua with only one root and in this case 
one has only field, canonical or ghost. When, however, more than one root is present ghosts 
appear in the spectrum. There are therefore two issues which should be clarified: what is 
the lifetime of such a vacuum before it gets destroyed by the presence of ghost and whether 
one can UV complete the model. 

As outlined above, on the one hand one might expect that considering higher order 
perturbations would make manifest ghosts instabilities that are not seen at the linear level. 
On the other hand, looking at (3.6), it is clear that the ghosts degree of freedom in our 
model is coupled not only through gravity with everything else, but also directly with x- 
The instability is due to the infinite phase space associated to the arbitrarily large negative 
energies of the ghost that result in a divergent decay rate for the vacuum. However, this 
instability might be not dramatic in an absolute sense. In fact, the phantom theory is meant 
to be just an effective description of certain regimes here. The validity of this description 
therefore depends on whether the ghost instability has the time to develop or not over the 
relevant cosmological time scale. 4 

We follow the analysis of [32] in order to give an estimate of the domain of validity of 
our specific model as an effective description with cutoff A. C utoff- I n [32] the allowed decay 
process involving phantom particles have been analyzed. We are particularly interested to 
the situation where there is just one ghost field, so as we are assuming in our case. Ordinary 
particles can decay into ghosts plus other ordinary particles, provided that the effective mass 
of the final ordinary particles is larger than the original one. When just one ghost is present, 
this can decay only if in the final states there are at least two ghosts and, at least, an 
ordinary particle. In order to compute the decay rate of our vacuum it is useful to rewrite 
the action involving the ghost field with the canonical normalization. We therefore perform 
the rescaling r/g — >■ r in the action (3.1) where, in the case under consideration, the index i 
denoting the different field run over one complex field r together with its complex conjugate. 
We then expand the scalar field Lagrangian in terms of the rescaled real field x an d ip, with 
r = 775 (x + i^) as in (3.6), obtaining 



xOx - if^Oip - mix 2 - V> 2 ) + 2nxV> 

- 3nG N (^( X 4 + ^) + ^^xV + xyx^X 2 ~ ^) 

with m = - ~ x + 6x -^° a nd n = v ^ x ■ To give a crude approximation on the lifetime of 
our vacuum, let us set to zero the quadratic term 2nxip- 5 Looking at the quartic part of the 
Lagrangian we see that several interactions involving phantoms are present. Following the 



4 See for instance [35] for an example where the phantom phase is just an intermediate phase in the 
cosmological evolution. Some unstable modes propagate out of this phase and some other not. 

5 This is actually a delicate point in our context. In fact setting n = is not just a simple choice of 
parameters. This would correspond to y = or x — 3Ho, implying Im(J) = 0. Therefore taking this specific 
limit do not correspond at all to a situation with two complex roots for the operator _F, but on the contrary 
corresponds to a situation with two real coincident roots. For simplicity let us however assume that we can 
neglect the quadratic term above. 
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rules above only a few decays are however allowed. In particular the kinematically permitted 
ones are ip — > ipipip, ip ~^ an d X ~^ XX^- The decay rate for a particle at rest for this 
kind of processes can be written generically as 

The prefactor is the mass (that is the same in all the processes here) of the particle that 
decay, M is the effective coupling constant of the quartic interaction corresponding to the 
process considered, in our case these are just tree level processes and the factor \A4\ 2 can 
be taken out of the integrals. The integrals over the momentum are up to a cutoff energy 
^■cutoff-, as explained above. The coupling constant of all the processes corresponding to 
quartic interaction in our Lagrangian is roughly just proportional to G7\r|a| 2 ~ \a \ 2 /Mp and 
all the phantom and the canonical field considered have the same mass squared, proportional 
to | J| ~ |a| 2 . On a purely dimensional analysis basis we can thus estimate the decay rate 
order of magnitude to be 

r (W\ 2 \ 2 K lutoff |«| 3 A 2 (7q x 
T72 —Z ~ 774 A cutoff • (7-3) 



Ml) rm M p cutoff ' 
If we consider the SFT origin of the class of model discussed, the order of magnitude of a is 



roughly of the same order as the string mass, which is the inverse of the string length 
and we can consider it to be related to Mp as a ~ zMp. Here we have introduced a new 
parameter z, which is just the ratio of the string mass to the Planck mass [36]. z may vary 
in a wide range since experimentally any mass above, say, 10 TeV is not excluded as a string 
scale. The natural cutoff coming from SFT is the string mass giving A cuto ff ~ zMp. The 
decay mean time for the processes considered above expressed in the unit of Hubble time , 
Hq 1 ~ 10 60 Af P , is of the order 

10 -60 

HoTdecay ~ ^ ■ (7.4) 

It is obvious from the latter that the lifetime of the vacuum can be large and even 
roughly comparable with the age of the Universe. The price for this is the lightweight string. 
z, which determines this, can be in the range 10~ 15 -j- 1, since we do not expect strings 
lighter than 10 TeV from experimental bounds and heavier than the Planck mass from the 
theoretical perspectives. 

Concerning the second question about the UV completion, the following idea might 
provide an answer in the context of SFT, even if the solution is most likely not unique. The 
point is that in the class of model we are considering here, in contrast to canonical scalar field 
models, different vacua can easily have different number and nature of degrees of freedom. 
One may therefore expect that, apart from the vacuum corresponding to the present model, 
there is another vacuum which has no ghosts and is generally well behaved. To illustrate this 
point we can consider the following action 

S = Jd 4 x Q0(D-m 2 )0 + ^(n)<^ , (7.5) 

where Q is a non-local operator and v describes some interaction potential. Now suppose that 
there are two minima with different second derivatives of the function v only. The kinetic 
operator would then look like 

T k = a-m 2 + Vk g{a), (7.6) 
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where Vf. = v"(G(0)(pk) and <fik are distinct vacuum solutions. A specific function of the form 
Q = «(□ — m 2 )(e^ D — 1) shows explicitly that in the vacuum with Vk = 1/a there is only one 
root for the resulting function T\. while any other value of generates infinitely many roots 
and therefore ghosts. Indeed, if v k = 1/a one has 

and this kinetic operator corresponds to one scalar degree of freedom. Even though this is 
not a genuine SFT derived action, it clearly illustrates how the model can be UV completed. 

The bottom line is that the parameters of the model can be such that an effect like the 
phantom divide crossing is observable and is extended over an appreciable lapse of time, but 
the eventual stabilization of the model is most likely achieved as a transition to a really stable 
vacuum and the discussion in the second half of this section describes how this can happen. 
Here it is important to mention that, as it is discussed in [33, 34], the ghost instability time 
scale can be finite only if the cutoff is Lorentz violating as in our consideration where we have 
used a preferred frame. This is not what one expects from the viewpoint of SFT. It is however 
not excluded that the false vacuum of SFT which contains ghosts could be Lorentz-violating. 
This question definitely requires better understanding. 

8 Conclusions 

In this paper we have studied a class of models which describe the late time evolution in the 
dynamics of scalar excitations in SFT. This class of models admits simple analytic solutions 
when minimally coupled to gravity. The main distinctive feature is the presence of non- 
locality arising from the presence of an operator ^(D) in the action. This non-locality can 
be translated in a (possibly infinite) set of local fields with masses squared Jj. These masses 
squared may be complex. The corresponding complex fields enter the Lagrangian in an 
unusual way, giving rise to a new structure. The main question we have addressed is the 
influence of such complex fields on the evolution of linear perturbations of the system. In the 
meantime, we have shown explicitly that these complex fields give a solution describing the 
crossing of the phantom divide. Therefore stability (or instability) under linear perturbations 
is essential to understand whether our model can (or cannot) be used as a viable model for 
the explanation of such a phenomenon. 

In order to have good background solutions we must require that the parameters Qj of 
the model satisfy Re(ctj) < 0. Then the solutions for the scalar fields do not blow up and the 
Hubble parameter tends to a constant. The most important result is the fact that, having 
such well behaved solutions, the perturbations related to complex oti eventually vanish for 
any values of the other parameters (the comoving wavenumber k, Im(aj), the initial data 
for the background solutions and perturbations). However, we notice that for some range of 
the parameters the perturbations do exceed the perturbative regime, making the question of 
validity of the linear approximation unclear. This range of parameters certainly deserves a 
more extensive analysis including the second order perturbations. 

From our analysis it turns out that in a general situation where real and complex Jj are 
involved, the perturbative spectrum is determined just by the real ones for a wide range of the 
parameters. Indeed, real values of Jj would generate ordinary real scalar fields. Cosmological 
perturbations with such scalar fields are well understood and are known to have growing or 
vanishing modes mainly depending on the values of the wavenumber k. This leads us to 
state that, since perturbations related to complex J; almost always vanish, the observable 
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spectrum of the model is determined solely by real Jj. Thus the crossing of the phantom 
divide, realized in the case of complex Jj, is perturbatively stable and perturbations do not 
grow even though the phantom phase is reached. 

The presence of ghost-like excitations however is alarming. We have computed an esti- 
mate for the vacuum decay rate due to ghosts in the context of our simple model. Admitting 
that the string mass is not necessarily of the order of Planck mass, one finds a quite large 
range of values for the lifetime. This can make the described effect of the phantom divide 
crossing observable. Here we want to stress again that, in SFT based models, the appearance 
of ghosts is most likely an indication of the fact that there exists another vacuum to which 
the model evolves. It is exactly the existence of a proper vacuum that is suggested as a way 
to make the model UV-complete. Even though there can be parameters such a that it is 
suitable for the present day observations, the model can have relevance for early stages of 
the Universe evolution, when effects like the phantom divide crossing were also possible. 

To summarize we conclude that models of type (3.1) provide interesting results from 
the point of view of cosmology. Apart from the phantom divide crossing, the appearance 
of non-gaussianity during the inflation is another fact which makes the class of models con- 
sidered in the present paper promising from the perspective of future cosmological studies. 
The results of this paper are generally applicable and one can state that perturbations in 
similar systems would not be triggered by new type of fields and quantum effects can be 
suppressed for a long time according to the vacuum decay rate estimate. Therefore, this 
new class of models is believed to provide new successful ways of resolving some of the cos- 
mological questions. Further aspects, like coupling to other cosmic fluids (CDM, radiation, 
etc), construction of other potentials, consideration of thermodynamical properties 6 as well 
as better understanding of the UV-completion mechanism certainly deserve a deeper study. 

Acknowledgments 

The authors are grateful to I.Ya. Aref'eva, F. Bezrukov, B. Craps, B. Dragovich, G. Dvali, 
V. Mukhanov and S.Yu. Vernov for useful comments and stimulating discussions. This work 
is supported in part by the Belgian Federal Science Policy Office through the Interuniversity 
Attraction Pole P6/11, and in part by the FWO-Vlaanderen through the project G. 0114. ION. 
A.K. is supported in part by RFBR grant 08-01-00798 and state contract of Russian Federal 
Agency for Science and Innovations 02.740.11.5057. 

References 

[1] K. Ohmori, "A review on tachyon condensation in open string field theories," 
arXiv:hep-th/0102085; 

I. Y. Arefeva, D. M. Belov, A. A. Giryavets, A. S. Koshelev and P. B. Medvedev, 
"Noncommutative field theories and (super)string field theories," arXiv:hep-th/0111208; 
W. Taylor, "Lectures on D-branes, tachyon condensation, and string field theory," 
arXiv:hep-th/0301094. 

[2] L. Brekke, P. G. O. Freund, M. Olson and E. Witten, "Nonarchimedean String Dynamics," 
Nucl. Phys. B 302, 365 (1988); 

P. H. Frampton and Y. Okada, "Effective Scalar Field Theory of p-adic String," Phys. Rev. D 
37, 3077 (1988); 

6 See, for instance, [37] about interesting and promising results concerning the thermodynamics of p-adic 
strings. 



- 16 - 



V. S. Vladimirov, I. V. Volovich and E. I. Zelenov, "p-adic analysis and mathematical physics," 
Ser. Sov. East Eur. Math. 1 (1994) 1; 

B. Dragovich, A. Y. Khrennikov, S. V. Kozyrev and I. V. Volovich, "p-Adic Mathematical 
Physics," Anal. Appl. 1 (2009) 1 [arXiv:0904.4205 [math-ph]]. 

[3] I. Y. Aref'eva, "Nonlocal string tachyon as a model for cosmological dark energy," AIP Conf. 
Proc. 826 (2006) 301 [arXiv:astro-ph/0410443]; 

I. Y. Aref'eva, "Stringy Model of Cosmological Dark Energy," AIP Conf. Proc. 957, 297 (2007) 
[arXiv:0710.3017 [hep-th]]. 

[4] N. Barnaby and J. M. Cline, "Large Nongaussianity from Nonlocal Inflation," JCAP 0707, 017 
(2007) [arXiv:0704.3426 [hep-th]]; 

N. Barnaby and J. M. Cline, "Predictions for Nongaussianity from Nonlocal Inflation," JCAP 
0806 (2008) 030 [arXiv:0802.3218 [hep-th]]. 

[5] G. Calcagni, "Cosmological tachyon from cubic string field theory," JHEP 0605, 012 (2006) 
[arXiv:hep-th/0512259]; 

G. Calcagni, M. Montobbio and G. Nardelli, "Route to nonlocal cosmology," Phys. Rev. D 76, 
126001 (2007) [arXiv:0705.3043 [hep-th]]; 

G. Calcagni and G. Nardelli, "Tachyon solutions in boundary and cubic string field theory," 
Phys. Rev. D 78, 126010 (2008) [arXiv:0708.0366 [hep-th]]; 

G. Calcagni, M. Montobbio and G. Nardelli, "Localization of nonlocal theories," Phys. Lett. B 
662, 285 (2008) [arXiv:0712.2237 [hep-th]]; 

G. Calcagni and G. Nardelli, "Nonlocal instantons and solitons in string models," Phys. Lett. 
B 669, 102 (2008) [arXiv:0802.4395 [hep-th]]. 

[6] N. Barnaby, T. Biswas and J. M. Cline, "p-adic inflation," JHEP 0704, 056 (2007) 
[arXiv:hep-th/0612230]. 

[7] T. Biswas, A. Mazumdar, W. Siegel, "Bouncing universes in string-inspired gravity," JCAP 
0603 (2006) 009. [hep-th/0508194]; 

T. Biswas, T. Koivisto, A. Mazumdar, "Resolution of the Cosmological Singularity in 
Non-local Higher Derivative Theories of Gravity," [arXiv:1005.0590 [hep-th]]. 

[8] A. S. Koshelev and S. Y. Vernov, "Cosmological perturbations in SFT inspired non-local scalar 
field models," [arXiv:0903.5176 [hep-th]]; 

A. S. Koshelev, S. Y. .Vernov, "Analysis of scalar perturbations in cosmological models with a 
non-local scalar field," [arXiv:1009.0746 [hep-th]]. 

[9] A. S. Koshelev, "Non-local SFT Tachyon and Cosmology," JHEP 0704 (2007) 029. 
[hep-th/0701103]. 

[10] N. Arkani-Hamed, S. Dimopoulos, G. Dvali et al., "Nonlocal modification of gravity and the 
cosmological constant problem," [hep-th/0209227]; 

A. O. Barvinsky, "Nonlocal action for long distance modifications of gravity theory," Phys. 
Lett. B572 (2003) 109-116. [hep-th/0304229]; 

S. Dcser, R. P. Woodard, "Nonlocal Cosmology," Phys. Rev. Lett. 99, 111301 (2007). 
[arXiv:0706.2151 [astro-ph]]; 

S. 'i. Nojiri, S. D. Odintsov, "Modified non-local-F(R) gravity as the key for the inflation and 
dark energy," Phys. Lett. B659 (2008) 821-826. [arXiv:0708.0924 [hep-th]]; 
S. Jhingan, S. Nojiri, S. D. Odintsov et al., "Phantom and non-phantom dark energy: The 
Cosmological relevance of non-locally corrected gravity," Phys. Lett. B663 (2008) 424-428. 
[arXiv:0803.2613 [hep-th]]; 

T. S. Koivisto, "Newtonian limit of nonlocal cosmology," Phys. Rev. D78 (2008) 123505. 
[arXiv:0807.3778 [gr-qc]]; 

C. Deffayet, R. P. Woodard, "Reconstructing the Distortion Function for Nonlocal Cosmology," 
JCAP 0908 (2009) 023. [arXiv:0904.0961 [gr-qc]]; 

S. Capozziello, E. Elizalde, S. 'i. Nojiri et al., "Accelerating cosmologies from non-local 



-17- 



higher-derivative gravity," Phys. Lett. B671 (2009) 193-198. [arXiv:0809.1535 [hep-th]]. 

[11] E. Komatsu et al. [WMAP Collaboration], "Five- Year Wilkinson Microwave Anisotropy Probe 
(WMAP) Observations :Cosmological Interpretation," Astrophys. J. Suppl. 180 (2009) 330 
[arXiv:0803.0547 [astro-ph]]. 

[12] M. Libanov, V. Rubakov, E. Papantonopoulos et al., "UV stable, Lorentz- violating dark energy 
with transient phantom era," JCAP 0708 (2007) 010. [arXiv:0704.1848 [hep-th]]. 

[13] N. Arkani-Hamcd, H. -C. Cheng, M. A. Luty et al, "Ghost condensation and a consistent 
infrared modification of gravity," JHEP 0405 (2004) 074. [hep-th/0312099]; 
N. Arkani-Hamcd, P. Creminelli, S. Mukohyama et al., "Ghost inflation," JCAP 0404 (2004) 
001. [hcp-th/0312100]; 

S. Mukohyama, "Accelerating Universe and Cosmological Perturbation in the Ghost 
Condensate," JCAP 0610 (2006) 011. [hep-th/0607181]; 

P. Creminelli, G. D'Amico, J. Norena et al., "The Effective Theory of Quintessence: the w<-l 
Side Unveiled," JCAP 0902 (2009) 018. [arXiv:0811.0827 [astro-ph]]. 

[14] C. Deffayet, O. Pujolas, I. Sawicki, A. Vikman, "Imperfect Dark Energy from Kinetic Gravity 
Braiding" JCAP 1010 (2010) 026. [c-Print: arXiv:1008.0048 [hep-th]] 

[15] E. A. Lim, I. Sawicki, A. Vikman, "Dust of Dark Energy," JCAP 1005 (2010) 012. 
[arXiv:1003.5751 [astro-ph.CO]]. 

[16] P. Creminelli, A. Nicolis, E. Trincherini, "Galilean Genesis: An Alternative to inflation," JCAP 
1011 (2010) 021. [arXiv:1007.0027 [hep-th]]. 

[17] V. Sahni, Y. Shtanov, "Brane world models of dark energy," JCAP 0311 (2003) 014. 
[astro-ph/0202346]; 

A. Lue, G. D. Starkman, "How a brane cosmological constant can trick us into thinking that 
W < -1," Phys. Rev. D70 (2004) 101501. [astro-ph/0408246]; 

A. S. Koshclev, T. N. Tomaras, "Towards a covariant model for cosmic self-acceleration," 
JHEP 0710 (2007) 012. [arXiv:0706.3393 [hep-th]]; 

M. R. Sctare, E. N. Saridakis, "Brancworld models with a non-minimally coupled phantom 
bulk field: A Simple way to obtain the -1-crossing at late times," JCAP 0903 (2009) 002. 
[arXiv:0811.4253 [hep-th]]; 

Y. Shtanov, V. Sahni, A. Shafieloo et al, JCAP 0904 (2009) 023. [arXiv:0901.3074 [gr-qc]]. 

[18] R. Emparan, J. Garriga, "Non-perturbative materialization of ghosts," JHEP 0603 (2006) 
028. [hcp-th/0512274]. 

[19] R. Kallosh, J. U. Kang, A. D. Lindc, V. Mukhanov, "The New ckpyrotic ghost," JCAP 0804 
(2008) 018. [arXiv:0712.2040 [hep-th]]. 

[20] A. Adams, N. Arkani-Hamed, S. Dubovsky et al., "Causality, analyticity and an IR obstruction 
to UV completion," JHEP 0610 (2006) 014. [hcp-th/0602178]. 

[21] G. Calcagni and G. Nardelli, "Non-local gravity and the diffusion equation," arXiv:1004.5144 
[hep-th]. 

[22] F. Galli, A. S. Koshelev, "Multi-scalar field cosmology from SFT: an exactly solvable 
approximation," Theor. Math. Phys. 164 (2010) . [arXiv:1010.1773 [hep-th]]. 

[23] A. S. Koshelev, "SFT non-locality in cosmology: solutions, perturbations and observational 
evidences," AIP Conf. Proc. 1241 (2010) 630-638. [arXiv:0912.5457 [hep-th]]. 

[24] I. Y. Aref'eva, L. V. Joukovskaya, A. S. Koshelev, "Time evolution in superstring field theory 
on nonBPS brane. 1. Rolling tachyon and energy momentum conservation," JHEP 0309 
(2003) 012. [hep-th/0301137]. 

[25] V. S. Vladimirov, Y. I. Volovich, "On the nonlinear dynamical equation in the p-adic string 
theory," Theor. Math. Phys. 138 (2004) 297-309. [math-ph/0306018]; 



- 18 - 



V. S. Vladimirov, "On the equation of the p-adic open string for the scalar tachyon field," 
[math-ph/0507018]. 

[26] L. Joukovskaya, "Rolling tachyon in nonlocal cosmology," AIP Conf. Proc. 957 (2007) 
325-328. [arXiv:0710.0404 [hep-th]]. 

[27] A. Anisimov, E. Babichcv, A. Vikman, "B-inflation," JCAP 0506 (2005) 006. 
[astro-ph/0504560]. 

[28] B. Feng, X. -L. Wang, X. -M. Zhang, "Dark energy constraints from the cosmic age and 
supernova," Phys. Lett. B607 (2005) 35-41. [astro-ph/0404224]; 

Y. -F. Cai, E. N. Saridakis, M. R. Setare et ai, "Quintom Cosmology: Theoretical implications 
and observations," Phys. Rcpt. 493 (2010) 1-60. [arXiv:0909.2776 [hep-th]]; 
H. Zhang, "Crossing the phantom divide," [arXiv:0909.3013 [astro-ph.CO]]. 

[29] V. F. Mukhanov, H. A. Fcldman and R. H. Brandenbcrgcr, "Theory of cosmological 

perturbations. Part 1. Classical perturbations. Part 2. Quantum theory of perturbations. Part 
3. Extensions," Phys. Rcpt. 215, 203 (1992). 

[30] J. M. Bardeen, "Gauge Invariant Cosmological Perturbations," Phys. Rev. D 22, 1882 (1980). 

[31] J. c. Hwang and H. Noh, "Cosmological perturbations with multiple fluids and fields," Class. 
Quant. Grav. 19, 527 (2002) [arXiv:astro-ph/0103244]. 

[32] S. M. Carroll, M. Hoffman, M. Trodden, "Can the dark energy equation - of - state parameter 
w be less than -1?," Phys. Rev. D68 (2003) 023509. [astro-ph/0301273]. 

[33] J. M. Clinc, S. Jcon, G. D. Moore, "The Phantom menaced: Constraints on low-energy 
effective ghosts," Phys. Rev. D70 (2004) 043543. [hcp-ph/0311312]. 

[34] R. P. Woodard, "Avoiding dark energy with 1/r modifications of gravity," Lcct. Notes Phys. 
720 (2007) 403-433. [astro-ph/0601672]. 

[35] Y. -F. Cai, T. Qiu, R. Brandcnberger, Y. -S. Piao, X. Zhang, "On Perturbations of Quintom 
Bounce," JCAP 0803 (2008) 013. [arXiv:0711.2187 [hep-th]]. 

[36] I. Ya. Aref'eva, A. S. Koshelev, "Cosmological Signature of Tachyon Condensation," JHEP 
0809 (2008) 068. [arXiv:0804.3570 [hep-th]]. 

[37] T. Biswas, J. A. R. Cembranos, J. I. Kapusta, "Thermodynamics and Cosmological Constant 
of Non-Local Field Theories from p-Adic Strings," JHEP 1010 (2010) 048. [arXiv:1005.0430 
[hep-th]]. 



- 19 - 



